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Twin boundaries are generic crystalline defects in noncentrosymmetric crystal structures. We 
study theoretically twin boundaries in time-reversal symmetric noncentrosymmetric superconductors 
that admit parity-mixed Cooper pairing. Twin boundaries support spin currents as a consequence of 
this parity mixing. If the singlet and triplet components of the superconducting order parameter are 
of comparable magnitude, the superconducting state breaks spontaneously the bulk time-reversal 
symmetry locally near the twin boundary. By self-consistently evaluating the Bogoliubov-de-Gennes 
equations and the gap functions we find two distinct phases: First, time-reversal symmetry breaking 
enhances the spin currents but does not lead to chiral supercurrents. A secondary phase transition 
then triggers a spin magnetization and an orbital supercurrent near the twin boundary. 



Initiated by the discovery of superconductivity in the 
noncentrosymmetric heavy Fermion compound CePtaSi, 
noncentrosymmetric superconductors have opened up 
new perspectives in the study of unconventional super- 
conductivity [IHS]. Due to the lack of inversion sym- 
metry, these systems feature antisymmetric spin-orbit 
coupling that breaks completely the SU(2) spin-rotation 
symmetry. As a consequence, the superconducting con- 
densate is a superposition of spin-singlet and spin-triplet 
Cooper pairs [9) [10] . The possibility of realizing different 
mixing ratios between singlet and triplet pairing is cru- 
cial to explore the topological properties of time-reversal 
symmetric (TRS) noncentrosymmetric superconductors. 
For example, quasi two-dimensional (2D) superconduc- 
tors of this kind have a Z 2 topological attribute when 
fully gapped (symmetry class Dili in the classification of 
Ref. [11]). The phases with dominant s-wave pairing {p- 
wave pairing) are topologically trivial (non-trivial) . They 
are separated by a gap-closing topological phase tran- 
sition [12] . Akin to the helical electronic edge states 
of a 2D Z 2 topological insulator, helical edge modes 
in the form of Andreev bound states transport a non- 
conserved spin current along the boundary of a non- 
centrosymmetric superconductor with dominant p-wave 
pairing [51I71IT5 15 j. Andreev bound states are a specific 
signature of unconventional Cooper pairing and directly 
manifest themselves in tunneling spectroscopy measure- 
ments and unique features of spin transport [5]. 

Very generically, the topology of phases of matter can 
be probed at defects, such as boundaries, lattice dislo- 
cations or vortices in a superconducting order param- 
eter |16j . In noncentrosymmetric materials, the crys- 
tal structure allows for another type of defect when 
two regions of space with the opposite inversion sym- 
metry breaking face each other in a single crystal. In 
fact, the formation of these so-called twin boundaries is 
unavoidable in the crystal growth processes of noncen- 



trosymmetric materials. A first step towards understand- 
ing the parity-mixed superconducting state at intrinsic 
twin boundaries has been undertaken in Ref. [6]. It re- 
vealed that TRS can be spontaneously broken at the twin 
boundary and that the twin boundary can host vortices 
with fractional vorticity. The latter offer an explanation 
for unusually slow flux line dynamics that has been ob- 
served experimentally in noncentrosymmetric supercon- 
ductors [T7]: Fractionalized vortices are trapped at the 
twin boundaries and form an efficient barrier against flux 
flow. 

In this Letter, we develop an extended phenomenol- 
ogy of the possible pairing states at the twin boundary 
by self-consistently evaluating the Bogoliubov-de Genncs 
(BdG) equations and the gap functions. We find two 
distinct phases of broken TRS at the twin boundary. In 
the first phase, the orbital degrees of freedom of the su- 
perconducting order parameter break TRS and enhance 
the helical spin currents along the twin boundary. Yet, 
contrary to naive expectations, no chiral current flows 
despite the broken TRS. This changes with a secondary 
transition to another TRS broken phase that features 
both a nonvanishing magnetization and an orbital super- 
current along the twin boundary. In our analysis, we 
clarify the relation between the bulk topological phase 
transition and the phase diagram of the states near the 
twin boundary, as well as the nature and spatial profile 
of the spin and charge supercurrents in each phase. 

We use a tight-binding model describing a 2D non- 
centrosymmetric superconductor with Rashba spin-orbit 
coupling that includes paring interactions based on spin 
fluctuations, in order to allow for the appearance of un- 
conventional pairing channels. We consider electrons 
hopping on a square lattice A of L x x L y sites r = 
(x, y) € A that is spanned by the orthogonal unit vec- 
tors a i , i — x,y. The Hamiltonian for the extended 
Hubbard model including spin-density interactions and 
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the Rashba spin-orbit coupling [ITS] reads 

H ■= - J2J2 C r+a lS ( t& ° - X r.a t ' &)s,s>Crs' 
v , i s 

+ \J Sr+at ' &r + &r,ai ' (^r+Oi x *$V)] 
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+ 17 X! 7 M ? M + V ^2 n r n r+a, , 
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where cj. s creates an electron with spin s = (f, \.) at 
site r g A. The antisymmetric spin-orbit coupling is 
a Rashba term of strength a r parametrized by \ r ai = 
ia r (z X Ui), i = x,y, where z is the unit vector normal 
to the plane of the lattice. The vector <r = (<r x , a v , a z ) 
denotes the three Pauli matrices and <r° the 2x2 
unit matrix. We define the electron density operator 
n Ti = n Ti ^ + n ri ^ and the spin density operator S r = 
Ess'4s^ss ,c ^s , ■ Besides the ordinary spin-isotropic 
Heisenberg exchange of strength J, the noncentrosym- 
metric crystal structure also allows for a Dzyalonshinsky- 
Moriya type spin-spin interaction of strength D r which 
is parametrized as D r a% = D r (z x a,), i — x,y. 

Let us illustrate the mean-field decoupling for a trans- 
lational invariant system (a r = a, D r = D, Vr), as- 
suming periodic boundary conditions in both the a x and 
a y directions. The Rashba-type spin-orbit interaction 
breaks the SU(2) spin-rotation symmetry and induces a 
splitting of the electron bands, with each band exhibit- 
ing a specific spin structure in momentum space. At the 
same time, it allows the gap function to be of mixed par- 
ity. Here, we focus on a scenario of mixed s- and a p-wave 
component for which the BdG mean-field Hamiltonian 
reads H BdG :— J^k c l^fc dGc fci wrtn tnc four-component 
notation Ck = (c^-j- , c^, cL fe -p cL fe j_) and k representing 
the 2D momentum in the Brillouin zone while 
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The kinetic part of the Hamiltonian is given by 



-^[cos^) + cos(/cj,)]ct° + <7fc ■ <t, 
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with Qk '■= 2a[xsm{k y ) — ysm(k x )]. We can decompose 
the superconducting gap function A fc = (ier y )(Aj^ + 
A^ ■ <x) in a scalar even parity spin singlet part A^ 

and a vector odd parity spin triplet part [HUB]. The 
assumption of s-wave pairing in the former and p-wave 
paring in the latter yields the momentum dependences 



A 



= (cos(k) + cos(fc,,)H 
= A p (sin(fc. y ), - sin(fc x ), 0) . 
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of the order parameters. Note that the p-wave paring 
state with the highest transition temperature is the one 
with A^ oc gk [Sj- This allows to simultaneously di- 
agonalize the gap function and the kinetic part of the 
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FIG. 1. (Color online) (a) Sketch of the intrinsic metallic in- 
terface in a noncentrosymmetric superconductor. If the pair- 
ing is dominantly of triplet type, a topologically protected 
Kramers pair of edge modes carrying a spin current is lo- 
calized at x = and x = L x . Two Kramers pairs of edge 
modes that are localized near the twin boundary hybridize 
and thereby loose their topological protection, (b) The spec- 
tral function of up-spin quasiparticles at the immediate left of 
the twin boundary , A^._ &0 , (E) shows the spin-polarization 
of these helical states. The zoom-in (c) reveals their hy- 
bridization. (A^ =50 k ^ (E) and A^ =51 k (E) are obtained by 
flipping the figure about k y = 0.) Parameters are J = 1.3, 
D = 1.75, V = 1.22, U = 0.82, such that A s /A p ~ 0.425 [15] . 



Hamiltonian 'Hkin.fc by going to the basis of A = ± helic- 
ity states that label the two sheets of the noninteracting 
band structure. The gap function on either sheet is then 
given by A A > = A^ c) + AA^ o) • g k /\gk\, A = ±, making 
the mixed parity character of the superconducting state 
explicit. 

We now turn to the electronic properties of the sys- 
tem with twin boundary. For that, we equip the Rashba 
spin-orbit coupling and the Dzyalonshinsky-Moriya inter- 
action in Hamiltonian (JlJ with the spatial dependencies 



asgn(x - L x /2), 



Dsga(x-L x /2), (3) 



respectively. This models a twin boundary located at 
x = L x /2 that separates regions A and B of the super- 
conductor, that have the opposite sign of the Rashba and 
Dzyalonshinsky-Moriya coupling [see Fig. [I] (a)]. Open 
and periodic boundary conditions are used in the a x and 
a y directions, respectively. In other words, the relative 
U(l) phase between A( c ) and A p has to change by tt 
across the twin boundary. The way the superconduct- 
ing condensate accommodates this phase twist decisively 
determines the physics at the twin boundary. 

If A' ) is the dominant component of the order param- 
eter, gapless helical edge states exist at the boundary 
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of a 2D topological superconductor, as dictated by the 
nontrivial Z2 topological index. Localized modes within 
the bulk spectral gap also exist at the twin boundary, 
though they would not be endowed with topological pro- 
tection, for the Z2 topological sector of the Hamiltonian 
is the same on either side. In order to identify the ex- 
istence and spin polarization of the localized modes at 
twin boundary, we calculated the spectral function 
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(E), 
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where G'\ is the Green's function at position x and 
momentum k y . Figures Qb) and (c) show the spectral 
function for up-spin quasiparticles on the immediate left 
of the twin boundary. We observe that the left- and 
rightgoing modes have opposite spin polarization at high 
energies and that this spin-momentum locking is lifted 
with the appearance of a hybridization gap near zero en- 
ergy. 

We are now prepared to report the central results of 
our analysis on the spontaneous breaking of TRS near 
the twin boundary. Time-reversal symmetry breaking is 
signaled by the following two quantities: (i) the relative 
U(l) phase 



(5) 



tp := argA^2 W4 - argA^ 3W4 

of the singlet component of the superconducting order 
parameter between regions A and B and (ii) the spin 
magnetization M cx — at the twin boundary. To 
understand the relevance of (i), we have to ask how the 
condensate can account for the 7r-shift in the relative U(l) 
phase between A^ and A p across the twin boundary. 
The two values of tp compatible with TRS are tp = and 
(p = 7T. In the former case, |A( C )| ^> |AW| and A p ac- 
counts for the relative phase shift of 7r with a sign change 
by going through zero at the twin boundary. Conversely, 
if |A( C )| < |A(°)|, A (c} will change sign and be nodal at 
the twin boundary, such that ip = n. In contrast, if |A( e )| 
and |A(°)| are of comparable magnitude, the system is 
frustrated and the cost in condensation energy for a node 
in either component at the twin boundary is high. In this 
case, the relative U(l) phase between A^ c ^ and A p will 
change continuously from to 7r across the twin bound- 
ary, thereby breaking the TRS locally. This, in turn, is 
signaled by ip 7^ 0, it, where either of the degenerate so- 
lutions (p and — tp is spontaneously chosen. The phase 
diagram in Fig. [2] shows that the local TRS breaking at 
the twin boundary as measured by tp (continuous lines) 
and the magnetization M (color code) result in two dif- 
ferent phase boundaries. We conclude that a secondary 
phase transition toward the phase labeled C with finite 
magnetization M arises inside the TRS breaking phase 
with ip ^ 0, 7r labeled B. 

For comparison, the dashed lines in Fig.[2ja) define the 
region in parameter space where the bulk system is gap- 
less due to nodes in the superconducting order parameter. 
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FIG. 2. (Color online) (a) Phase diagram of superconducting 
and magnetic order at the twin boundary as a function of 
on-site repulsion U for J = 1.3, D = 1.8, V = 1.2 [!§]. Con- 
tinues lines separate phases of (A) TRS singlet dominated 
superconductivity, (B,C) TRS breaking phase where the tp 
deviates from the values and n near the twin boundary, and 
(D) TRS triplet dominated superconductivity. The color scale 
measures the spin magnetization near the twin boundary that 
becomes nonzero in phase (C) via a secondary phase transi- 
tion in the phase with broken TRS. Dashed lines encircle the 
region in parameter space where the bulk superconducting 
gap is nodal, marking the topological phase transition from 
the trivial (left) to the nontrivial (right) Z 2 topological super- 
conductor, (b) Spin-magentization M and phase tp along the 
dotted line in (a). 



Left (right) of this region, the bulk is a topologically triv- 
ial (non-trivial) superconductor. The topological phase 
transition acquired a finite width due to the term A s in 
the Hamiltonian. Inside the gapless region, the order pa- 
rameter has several point nodes on the Fermi surfaces. In 
the edge Brillouin zone, flat edge bands stretch between 
the projections of pairs of these nodal points, much like 
the surface flatbands in graphene zigzag edges. We make 
the following two observations that relate to these topo- 
logical features. First, the bulk topological phase tran- 
sition happens fully inside the parameter range in which 
TRS is spontaneously broken at the twin boundary. Sec- 
ond, the magnetic order at the twin boundary is strongly 
enhanced in the bulk gapless region. This can be at- 
tributed to Stoner-like magnetism of the edge flatbands 
at the twin boundary. 

Finally, let us study how the spontaneous TRS break- 
ing manifests itself in the (spin) Hall response at the twin 
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FIG. 3. (Color online) (a) Spin current in the y direction as 
a function of x. The solid lines and the dashed lines are the 
spin currents with spin polarization in the x and y direction, 
respectively. Parameters are J = 1.1, D = 2.0, V = 1.5, U = 
0.87, such that A s /A p = 0.41 1!) . (b) The x-dependence of 
the orbital supercurrent (red line) and of the magnetization 
(blue line) near the twin boundary in the TRS broken phase. 
Parameters are J = 1.3, D = 2.1, V = 1.2, U = 0.87, such 
that A s /A p = 0.5 p2]. 



boundary. For that, we define the spin current of polar- 
ization i = x, y, z that runs in the y-direction 



rspin— i 

Jy 



(6) 



proximity to a topological phase transition can combine 
in a nontrivial way to create a phase of spontaneously 
broken TRS at the interface. The twin boundary either 
hosts a spin current or carries a chiral supercurrent along 
with a spontaneous spin magnetization. The TRS break- 
ing phase is intimately connected with the bulk topolog- 
ical phase transition. 

In closing, we note that the mechanism for TRS break- 
ing that we discussed in this work is not limited to twin 
boundaries in noncentrosymmetric superconductors. It 
can also be of relevance to tunable devices with interface 
superconductivity such as SrTiC>3 /LaAlC>3 , if the Rashba 
spin-orbit coupling is not uniform in space [21 . A similar 
phenomenology might apply to other ordering phenom- 
ena, such as the spontaneous generation of the quantum 
spin Hall effect in graphene-like materials [22] . In this 
case, TRS may be broken at the boundary between re- 
gions of opposite spin-Hall conductivity, spontaneously 
generating a charge Hall effect. 

This work was supported in part by the Swiss National 
Science Foundation. E. A. was supported by a Grant-in- 
Aid from JSPS. 



where the trace is taken over all states below zero en- 
ergy. We show in Fig. [3^a) the spin currents Jy plu ~ x 
and Jy Pin ~ z as a function of position x inside the phase 
B of phase diagram [2]ja). Each of them has opposite 
signs on either side of the twin boundary. The compo- 
nent Jy pm ~ v vanishes in the bulk and is much smaller 



than J^p 
component Jy P 



at the twin boundary (not shown). The 
corresponds to the usual spin current 
carried by an Andreev bound states and is similar at the 
edge of the sample and at the twin boundary. The com- 
ponent Jy pin ~ x is finite also in the bulk and increases 
at the boundary. However, its bulk contribution should 
not be interpreted as a physically measurable spin cur- 
rent [5D] . The form of the boundary contribution of the 
spin current is governed by the length scale that the spin- 
orbit coupling a introduces, namely the inverse differ- 
ence between the Fermi momenta of the spin-orbit split 
Fermi surfaces 1 / 1 — fcp_ | . If this length scale becomes 
smaller than the coherence length that governs the decay 
of the spin current towards the bulk, oscillations appear 
in Jy pin ~ x (not shown in Fig. [3]). The phase C of phase 
diagram]^ a) is characterized by the finite magnetization 
M shown in Fig. [3^b). This magnetization is accom- 
panied by a orbital supercurrent that runs in opposite 
directions on the immediate left and right of the twin 
boundary [Fig. (3^b)]. 

In summary, we studied the superconducting interface 
states between twin domains in a mixed-parity super- 
conductor by evaluating the BdG equations and the gap 
functions self-consistently. We find that the breaking of 
inversion symmetry, the superconducting order, and the 
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